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Abstract 

We report a class of integrable one-dimensional interacting electronic systems 
with off-diagonal disorder. For these systems, the disorder can be "gauged 
away," and the spectrum can be mapped completely onto the spectrum of the 
ordered problem, which can be solved by Bethe ansatz or by bosonization. We 
study the magnetic properties of the persistent currents in mesoscopic rings 
in the case in which the ordered system is a Luttinger liquid. The system can 
be paramagnetic or diamagnetic, depending on the amount of the disorder 
and the number of fermions in the system. 
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It is usually believed that any amount of disorder in one- dimensional systems leads to 
electron localization. This belief is based mainly on the work of Anderson [I] and Mott [0, 
which demonstrates that in a one dimensional, non-interacting electronic system, disorder 
leads to destructive interference and as a consequence localization of the electronic wave 
functions. Recently, a class of 1-dimensional models has been discovered in which at certain 
energies, localization is suppressed as a result of a subtle resonant effect. The random dimer 
model ||, used for instance to explain the high conductivity of doped polyaniline, ||, is 
the archtypeal model here. When interactions are included in the presence of disorder, 
even more complex behaviour can arise. In some cases, such as v = 2/3 edge states in the 
fractional quantum Hall effect, the disorder appears to act as an irrelevant operator ||. In 
addition, it has been shown at the level of mean-field theory that interactions might destroy 
resonant scattering in 1-dimensional systems ||. This issue remains controversial as only a 
mean-field answer has been obtained. 

In this paper we present an integrable, interacting model with disorder. We consider 
a particle hopping on a one-dimensional chain with a hopping matrix element that has a 
random phase. That the matrix element is complex means that time- reversal symmetry 
has been broken. Such a situation can arise, for instance, in the presence of strong spin 
orbit coupling [[?]]. A similar type of problem was studied in the late seventies by Mattis 
|| and by Luttinger || in the context of the Ising model. Apart from its intrinsic interest 
as an integrable system, our model is relevant to understanding the problem of persistent 
currents in mesoscopic rings. As we will show, the zero temperature magnetic properties of 
the disordered integrable system are quite remarkable: the system can be paramagnetic or 
diamagnetic, depending on both the number of fermions (an effect that has been studied 
previously |10[) and the amount of disorder in the system. 

We will consider the problem of interacting electrons with charge e moving in a ring 
(periodic boundary conditions) with N a sites which is subject to an external magnetic flux 
$. The Hamiltonian of interest is 
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N a ^ 

H = - J2 hi+ie 7 ^ (4,aCi+i,a + cl +ha c iitT ) + V[n ha ] (1) 

i=l,<7 

where c iiCr and cj jCr are the annihilation and creation operators for an electron in the i th site 
with spin projection a (a =| or |), V[n ii(7 ] describes the interaction between the electrons 
(which is assumed to depend only on the electron density, i.e., number of electrons on a 
site (rij )<7 = c i a Ci :a )), Uj+i is the hopping matrix element (which will be the source of the 
disorder), and $ = 7 is the flux quantum. This Hamiltonian is quite general in terms 
of the electron-electron interaction; for instance, it can describe Hubbard models involving 
only an on-site interaction, V = V\ = U^ni^n^i, extended Hubbard models having also 
nearest-neighbor interactions, V = V\ + V J2i,a n i,a n i+i,a, etc.. The disorder in the model 
arises from the hopping matrix element, which is chosen to be, 

hi+i = te ik (2) 
where t is a constant and Si is a random phase on the site which can vary between — n < 

Si < 7T. 

Becasue the disorder appears as just a phase in the Hamiltonian, it can be (almost !) 
completely "gauged away" by an unitary transformation, 

N a N a n-1 , 27r$ 
3=1 n=l,a 



U = exp { -i £ £ E ( Si + j^-J c)^ } (3) 



which transforms the fermion operator to 

d j>a = Uc 3 , a U ] = exp ji ^£ 5 i + U ~ ^jtf^J I c i>°- ^ 

Under this unitary transformation it is easy to see that the new Hamiltonian can be written 

as, 

H' = UHW = -tJ2 (4, ff di+i,a + 4+i,A<0 + V[v], (5) 

i=l,cr 

since Urij j(7 W = n^ c . Eqn. (5) describes the underlying ordered problem, simply reexpressed 
in terms of the dj >a ; thus we conclude that (modulo boundary conditions, as we discuss below) 



the integrability of the disordered problem depends only on the integrability of the ordered 
problem (5). Here we will consider cases in which Eqn. (5) can be solved by the Bethe 



ansatz [O or by bosonization [12 



Because we are treating a ring, we choose periodic boundary conditions, that is, Cj tU = 
Cj + N a)CT . In terms of the new operators, this requires 

d j+Na ,a = exp jz (A + 2vr-J^ j d jt (6) 

where A = Y^=x 5 t . Note that A can be always written as 

A = 27m + 5 (7) 

where n is an integer and — it < 5 < it. Observe that the effect of the disorder is to create 
an effective flux, $ e //, through the ring. This flux is given by 

$ ef f = $ + n$ + ^. (8) 

Thus, disorder in this problem is effectively equivalent to an increase in the magnetic field 
applied in the ring. Without loss of generality, we will assume that n = in what follows. 

It is interesting to look at the wavefunction of the original disordered problem, >, in 
terms of the wavefunction of the underlying ordered problem, >, which is an eigenstate 
of the Hamiltonian (5). Under the unitary transformation (3) we have, 

|*>=t/ t |* / >. (9) 
The general form for the wavefunction for the ordered system is 

JV 

Ul,j2, • • • , 3n) 4i,<Ti d k«^ ' ' ' 4n,*n i° > ( 10 ) 

where N is the number of electrons in the lattice, |0 > is the vacuum state of the problem. 

Let us consider first the case where (5) can be solved by the Bethe ansatz with the 
twisted boundary conditions (6). In this case the wave amplitude ip can be written as [11| 



^ 1 , ff2 ,.,. iV (ji,j2,---,J J v)=E^ 1 ^,---^(2|^)e iEi -^ : ' i (11) 

v 



for a specific ordering {1 < jox < 3Q2 < ■■■ < Jqn < N a } of the particles. V and Q are 
permutations of the coordinates ji and momenta k\ of the electrons (we measure momentum 
in units of the inverse lattice constant) and 

A^^MV) = (-l) r '^. Sl ,. S2 ,.,. ajV (Q|^) (12) 

are the amplitudes obtained by the Bethe ansatz equations for the ordered case (the 
factor ifp is or 1 if the permutation is even or odd, respectively). 

The effect of disorder in the wave function can be obtained from (9) by introducing the 
identity operator between the operators in (10) (UW = 1). It is straightforward to show 
that the wavefunction of the disordered system is simply written as (<3> = 0), 



N 



l*>=IIE ^ 1 ,a 2 ,.., CTJV (jl,j2, • • • ,jjv)cJ liffl ct 2iff2 • • -ct^JO > (13) 
2=1 JIM 



where 



(ji,j2, • • ■ , jv)exp ^7r^iV^J ( 14 ) 

where N t — N — m for j Qm < / < jg( m +i) — 1 in the ordering defined above. Thus, we 
conclude that the effect of this particular kind of disorder is to introduce a Berry phase in 
the wavefunction ||13| . Therefore, the disorder will introduce fluctuations in the shape of the 
wavefunction. For example, for the non-interacting problem, instead of an extended Bloch 
wave, we would find an extended but strongly oscillating solution. Of course the solution 
of the problem now depends on the specific form of the Hamiltonian which is solved by the 
Bethe ansatz with the twisted boundary conditions. 

As an important example of the relevance this model, let us consider the transport 



properties of this system at low energies. Using a well-known result due to Kohn |14|, we 
write the charge stiffness of this system as 



flc = yUrl (15) 



4>=0 

where Eq(^>) is the ground state energy of the system in the presence of the magnetic flux. 



Eq(&) can be calculated from the Bethe ansatz solution |15|Jl6l| . To explore the implications 



for transport in detail, let us take the continuum limit of the Hamiltonian (5), which (for 
typical band fillings) is a Luttinger model. For simplicity, let us consider the specific case of 
spinless fermions which interact via 

V[m] = u£(n-l) U + i-l). (16) 



i=i 



This problem has been solved for an external flux ]I0| , using the bosonization technique 



with twisted boundary conditions. We can apply the results of |H| directly, substituting 



for the external flux the effective flux due to disorder, as defined in (8). With some trivial 
modifications of the formula for D c in ref. [10] one finds, at zero temperature and for an odd 



number of electrons, 



^ = ^°£(-ir +lCOsM 



„=i n 



D° c In (2 cos 5 -\ (17) 



and for an even number of electrons 



n=l U 



= D° C In (2(1 -cos 8)) (18) 

where D® = ^ and v* F is the velocity of propagation of the fermions renormalized by the 
interactions For an odd number of electrons, the system is diamagnetic if \8\ < ^ and 
is paramagnetic if ^ < \8\ < it. For an even number of electrons, the system is diamagnetic 
for I < \S\ < 71 and paramagnetic for \8\ < |. Thus, the magnetic response of the system is 
very sensitive to both the "parity" (N) and the amount of disorder (5) in the system. 

The above argument is valid for a fixed value of 5. Suppose, however, that we have an 
ensemble of rings with different values of the disorder. We asume now that 5 is randomly 
distributed between —a and a where a < it, that is, the probability of finding an amount 
of disorder between 8 and 8 + d8 is given by 

P(8)d8= — Q(a- \8\)d8 (19) 
2a 



where 6 is the usual Heavyside step function. In this case, from (8) we would find that the 
average effective flux in the ring equals the applied flux, that is, $ e // = $. It is easy to 
show by averaging (17) and (18) with (19) that, for an odd number of electrons one finds, 

D c = D° c - J2 { -^^smna (20) 

a n=l U 

and the response is always diamagnetic and for an even number of electrons, 

1 00 1 

D c = -/J C °-]T— sinw* (21) 

and the response is always paramagnetic. This behavior resembles the problem without 
disorder. However, it is easy to see that for a = ir the charge stiffness vanishes in both 
cases. This is due to the random distribution of diamagnetic and paramagnetic rings which 
gives rise to an effect that averages to zero. This is clear in (17) and (18) since the stiffness 
is periodic in 5. 

In summary, we have presented in this paper a class of integrable disordered one- 
dimensional (periodic) systems in which the disorder can be "gauged away" . In such systems, 
the integrability only depends on the integrability of the underlying ordered problem. For 
problems solvable by the Bethe ansatz, the wavefunction acquires a Berry phase due to the 
motion of the electrons in the disordered system. We show that the disorder acts as an ef- 
fective flux through the ring and that, for a fixed total amount of disorder, beyond the usual 
sensitivity of the persistent current to the number of electrons in the system, the magnetic 
response depends also on the amount of disorder in the system. 
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